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We show that the quantum fluctuations associated with 
the Bogoliubov quasiparticle vacuum can be strongly concen- 
trated inside dark solitons in a trapped Bose Einstein conden- 
sate. We identify a finite number of anomalous modes that 
are responsible for such quantum phenomena. The fluctua- 
tions in these anomalous modes correspond to the 'zero-point' 
oscillations in coupled dark solitons. 



PACS numbers: 03.75.Fi, 05.30.Jp 



Dark solitons in quasi one-dimensional Bose-Einstein 
condensates (BEC) are coherent structures created from 
macroscopic excitations. Recent experiments have real- 
ized such dark solitons by using the phase imprinting 
technique , and several theoretical studies have ana- 
lyzed the stability and dynamical properties , as well 
as novel soliton-pair solutions in two-component Bose 
condensates |^||. Like other familiar topological states 
in BEC, dark solitons are described by macroscopic wave- 
functions. Even though atoms are excited, most of them 
still share the same wavefunction. Such a quantum de- 
generacy is the basis of treating the many-particle system 
as a single coherent matter wave. This raises a fundamen- 
tal issue about the robustness of quantum degeneracy as 
a BEC is excited to higher energy states that typically 
contain more than one solitons. The presence of multi- 
ple dark solitons and their mutual interaction would in- 
troduce extra quantum fluctuations affecting the macro- 
scopic coherence. In this paper we investigate the mech- 
anism of such quantum decoherence effects. 

We approach the problem by studying the quantum 
fluctuations associated with a chain of dark solitons in 
a trapped condensate. Specifically, we are interested in 
the minimal fluctuations represented by the quasiparti- 
cle vacuum in the Bogoliubov theory. In the case of 
ground state condensates, the corresponding quasipar- 
ticle vacuum leads to a cloud of incoherent atoms (usu- 
ally called condensate depletion) not described by the 
Gross-Pitaevskii mean field equation Such incoher- 
ent atoms are typically a small fraction of the condensate, 
their effects can be ignored in most situations involving 
dilute condensates. However, we find that the quasiparti- 
cle vacuum associated with the dark solitons has distinc- 
tive features accessible for experimental observations. In 
this paper we show that incoherent atoms can be strongly 



concentrated at the notches, and in addition we discover 
that this phenomenon is mainly the quantum fluctuations 
in a finite number of collective modes. The particle- like 
behavior of dark solitons is the key to interpret our find- 
ings. We shall see how dark solitons under confinement 
couple together as particles to execute collective oscilla- 
tions. Quantum uncertainties in both the positions and 
momenta impose a fundamental limit of the degree of 
'darkness' inside the dark solitons. 

To begin we consider a Bose-Einstein condensate of 
weakly interacting atoms confined in a harmonic poten- 
tial. A quasi one-dimensional condensate can be achieved 
in a trap in which the transverse motion is tightly con- 
fined (frozen) and only the longitudinal motion remains 
active. Given that m is the atomic mass and lot is the 
trap frequency, the second quantized system is modeled 
by the Hamiltonian (in units of %ujt), 
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Here the spatial coordinate x is in units of (h/muix) 
and f\ is the interaction strength between atoms averaged 
over the transverse area | pT| . As usual, the field operator 
^(x, i) can be decomposed into a dominant coherent part 
and a fluctuation part, i.e., \fr = y/~N $>+(/), where N is the 
number of condensate atoms. The mean field wavefunc- 
tion $ = $„e~ v "* is the n-th excited state of the BEC 
satisfying the nonlinear Schrodinger eigen-equation, 
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where fi n is the eigen-energy and g = Nfj. The quan- 
tum number n is the number of nodes appearing in . 
Therefore <&„ represents a static n-dark soliton solution. 

In Fig. 1, we show the typical shapes of $ n obtained 
from exact numerical solutions for the case g = 200. Al- 
though there are no simple analytic expressions for 
we find that variational methods are useful to approxi- 
mate these wavefunctions. We construct the trial wave- 
function as a product of single dark soliton solutions em- 
bedded in a smooth background with the ground state 
profile, i.e., 



$ n {x) ~ a n $ (x) JJtanh 
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Here a n is a normalization constant, qj is the position 
of the j-th dark soliton, and l^ 1 = y 1 / 2 <i>o(0) is the 
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inverse of the healing length (neglecting the slight po- 
sition dependence near the trap center). The ground 
state wavefunction <&o(x) can be further approximated by 

g~ 1//2 (p>o — x 2 /2) 1 ^ 2 in the Thomas-Fermi limit g ^> 1. 
We treat qj as variation parameters such that the mean 
field energy 
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is stationary. With the trial wavefunction given in Eq. 
(3), the qj obtained have about 1% error compared with 
the exact numerical results in Fig. 2. Such a good agree- 
ment is understood from the fact that dark solitons are 
sufficiently separated from each other by several healing 
lengths Iq, which allows the product form in (3) to be a 
close approximation in the Thomas- Fermi limit. Indeed, 
we find that the separation between the solitons is be- 
tween 3/o to 6Zo for a wide range of interaction strengths 
10 2 < g < 10 3 . We remark that dark solitons repel each 
other p2| . The mechanical equilibrium is attained by the 
inhomogeneous condensate background that provides a 
confinement 'force' holding the solitons together j5j. 

Having described the n-dark soliton mean-field solu- 
tion, we now proceed to discuss the fluctuation part 
4> n = 3' — VN<& n e~ itJ,nt . Following the standard Bogoli- 
ubov method [ft3|], we write </>„ as, 
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Here Cj and Cj' l> are quasi-particle creation and annihi- 
lation operators associated with the j-th collective mode. 
The mode functions u^ n \x), Vj n \x) and frequencies u>j 
are determined by 
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where £ = | 
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restriction: J dx([ 



v 2 + 4g |$„| - 2(i n j , under the 

4™^l 2 — 1^ | 2 ) = f - Using the expan- 
sion given in (5), the Hamiltonian (1) is approximated 
by a diagonalized form: H m E n + lu^c^ 1 ^ Cj 1 ^, where 

cubic and higher terms of <p n are ignored. 

Our task is to determine the incoherent part of atom 
density 



p n {x,t) = ($l(x,t)4> n (x,t) 
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associated with the vacuum state in which there are no 
quasiparticles in all collective modes, i.e., ( c^-Cj ) = 
0. It is not difficult to show that among all physi- 
cal states, the vacuum gives the minimal time-averaged 
p n (x,t). With the help of Eq. (5), p n (x,t) arising 



from the vacuum fluctuations is given by: p n (x) = 
J2j \ v j ( x )\ 2 - An integration of p n (x) over x gives the 
total number of incoherent atoms SN, i.e., SN = J2j 

with SNj = fdx\vi n) (x)\ 2 . 

We solve Eqs. (6) and (7) numerically for the first 
thousand collective modes in order to study the distribu- 
tion of incoherent atoms. Although SN involves all col- 
lective modes, we find that only the low-frequency modes 
are significant. This is shown in Fig. 2 in which SNj as 
a function of u>j is plotted. For the cases with soliton 
numbers up to four, we see that SNj diminishes quickly 
as the frequency increases. If we superimpose the four 
cases in Fig. 2 together, the four tails associated with 
the high frequency part of SNj all roughly lie on the 
same curve. This implies that the presence of more dark 
solitons essentially do not cause extra fluctuations in the 
high frequency modes. 

From now on we will focus on low frequency modes. 
One prominent feature we discover is that the number 
of modes with negative frequencies equals the number of 
dark solitons. These particular modes are shown in Fig. 
2 as black filled circles. In the Bogoliubov theory, collec- 
tive modes with negative frequencies and positive norm 
(i.e., J dx(\u^\ 2 — \ 2 ) — 1) are known as 'anoma- 
lous modes'. Previous studies on a single dark soliton 0] 
and a single vortex |H have identified the particle-like 
motion associated with anomalous modes. For example, 
in the n — 1 dark soliton case, the anomalous mode de- 
scribes the small amplitude oscillations of the soliton at 
a frequency —lut/V2 in the Thomas- Fermi limit ^,^,|). 
This agrees with the exact anomalous mode frequency 
shown in Fig. 2a within 1.3 percent. By increasing the 
interaction strength g, we find that the numerical fre- 
quency approaches —ujt/V2- 

In Fig. 3 the incoherent atom density p n (x) is plot- 
ted. We see that p n (x) (solid lines) are strongly concen- 
trated at the positions of the dark solitons. In each fig- 
ure the contribution from the anomalous modes is shown 
by the dash line, which is the partial sum in Eq. (9), 
counting only the finite number of modes with u)j < 0. 
By comparing the curves, we conclude that incoherent 
atoms near the dark solitons mainly originate from the 
vacuum fluctuations in the anomalous modes. In fact, if 
we remove the background effect by defining the quantity 
A(x) = p n (x) — po{x) as a measure of fluctuations due to 
the presence of dark solitons, we find a more striking ev- 
idence. This is shown in the inset of Fig. 3 in which the 
solid lines (A(x)) and the dash lines (from the anomalous 
modes) are almost indistinguishable from each other. 

Having elucidated the significance of anomalous 
modes, we now make a connection to the particle-like 
motion of the dark solitons. The fact that the number 
of anomalous modes and the number of dark solitons are 
equal (see Fig. 2) provides an important clue. We can 
imagine the n mutually repulsive dark solitons, each be- 
ing pushed toward the trap center by the inhomogeneous 
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background, form a one-dimensional 'lattice' at equilib- 
rium. As a result, there are n normal modes of oscilla- 
tions when the equilibrium system is slightly disturbed. 
These modes correspond to the anomalous modes that 
are well localized around the dark solitons. For example, 
the frequency ~ —lut/V^ appears in all four cases in Fig. 
2. We find that this mode corresponds to the in-phase 
mode that the n solitons move together without changing 
their relative separations. 

To elaborate the physical picture quantitatively, we 
present a method to estimate the effective masses and 
frequencies of the coupled oscillations of solitons. For 
dcfmitcncss, let us examine the n — 2 case, and consider 
the ansatz wavefunction of the dark soliton pair, 



$ 2 (z, t) = a<S> (x)ipi{x,t)Tp 2 {x,t). 



(9) 





Here a is a normalization constant and ipj(x,t) is the 
exact solution of a dark soliton at the position qj moving 
at a speed qj (j — 1,2) in a homogeneous background, 
i.e., 



V s 
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where v s = g 1 / 2 $o(0) is the local speed of sound for the 
ground state at the origin . The construction of $2(2:, t) 
generalizes the trial wavefunction (3) by including the 
speed qj as variational parameters. Our strategy is to 
determine how the energy of the system changes if qj and 
qj slightly deviate from the stationary point: qj = and 
qi = — q 2 = Q. This can be done by inserting the ansatz 
(9) into Eq. (4) and keep only the quadratic terms. The 
energy change can be cast in the form of two harmonic 
oscillators, 
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Here we have defined the generalized coordinates for the 
in-phase mode (+) and the out-of-phase mode (— ) by: 
Q± = ($<li ±8q 2 )/V2 with Sqi = q% — Q and Sq 2 = q 2 + Q- 

We find that ui± obtained in Eq. (11) match quite well 
with the exact numerical values of the anomalous mode 
frequencies, noticing that ansatz (9) is not an exact so- 
lution. For example, uj + and w_ attain 93% and 96% 
accuracy respectively for the case with g = 200. The ef- 
fective mass M for each soliton turns out to be negative, 
which is a basic property of dark solitons discussed in 
literature Q. When g ^> 1 is sufficiently large, we find 
M m -4ZV* (0) 2 Zo ~ -3.2Ng- 2 / 3 (in units of the atom 
mass m). Since g — Nrj is proportional to N, the abso- 
lute value of the soliton mass M increases as jV 1 / 3 as the 
number of particles increases. 

Treating the positions qj and momenta pj = Mqj as 
dynamical variables, we construct a two-particle Hamil- 
tonian based on the energy (11). The quantized form of 



the Hamiltonian obtained corresponds to the part of the 
Bogoliubov Hamiltonian containing the two anomalous 
modes. Therefore the creation and annihilation opera- 
tors associated with the anomalous modes can be explic- 
itly related to the position and momentum operators <jj 
and pj. In this way the vacuum state of quasiparticles 
in the anomalous modes are now recognized as the zero- 
point oscillations of the coupled dark solitons jn| . 

For the case n = 2 discussed here, the vacuum state 
gives the position fluctuations (Ag 2 .) = 1/2Muj±. Since 
the width of each dark solitons is about Iq, the fluctu- 
ations become significant near the notches if the ratio 
(Aqj) /Iq (j = 1,2) is greater or comparable to one. We 
find that this occurs if N 1 / & fj 2 /' i > 1, according to the 
expression of M obtained above and the fact that uj± are 
order one. We can estimate the number density of inco- 
herent atoms at the notches. For example, by counting 
only the contributions from the two anomalous modes, 
we have Jlql 



l 2 {x = ±Q) 



V2 



■ P+±P- q+±Q- 
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Therefore the density \ft>\{Q)(j> 2 {Q)j is approximately 

3 4 /3 2 -i/3£jVM-y /3 (w^+wZ 1 ) /32, which is an in- 
creasing function of g. Such an analytical estimation 
matches the exact numerical results in Fig. 3b with about 
92% accuracy. If a higher accuracy is desired, the ansatz 
(9) should be generalized. The exact solution of a dark- 
soliton pair by Blow and Doran ]l2[| in a homogeneous 
background should provide insights for making correc- 
tions. 

To conclude, we have investigated the fundamental 
fluctuations in dark solitons arising from their Bogoli- 
ubov quasiparticle vacuum. These quantum fluctuations 
are revealed in the form of incoherent atoms strongly con- 
centrated near the notches. We interpret such a quan- 
tum effect as a consequence of zero-point oscillations of 
the coupled dark solitons. By constructing a particle-like 
dynamical model with solitons' positions and velocities as 
(finite numbers) degrees of freedom, quantitative predic- 
tions can be made approximately. In particular we find 
the regime iV 1 / 6 ?^ 2 / 3 > 1 where the position uncertainty 
becomes comparable with the soliton's width. Therefore 
blurred dark solitons would occur if the particle number 
is high or the transverse area of the condensate is suffi- 
ciently small. 
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FIG. 2. The number of incoherent atoms associated with 
each collective excitation mode with the frequency uij. The 
four plots correspond to the first four dark soliton states in 
Fig. 1 with the same interaction strength g — 200. The 
negative frequencies modes are anomalous modes plotted in 
black solid circles. The frequencies ujj are in units of lot- 
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FIG. 3. Spatial density profiles of incoherent atoms (solid 
lines): (a) n = 1, (b) n = 2, (c) n = 3 and (d) n = 4 
with the same parameters used in Fig. 1. The dash line 
in each figure indicates the contribution from the anomalous 
modes only. The solid line in each inset shows the quantity 
A = p n {x) — po(x) as a function of a;. In each inset we have 
also plotted the contribution from the anomalous modes in a 
dash line (that appears almost exactly the same as A). 



FIG. 1. The first four stationary dark soliton states ob- 
tained from the numerical solutions of Eq. (2) with g = 200: 
(a) n = 1, (b) n — 2, (c) n = 3 and (d) n = 4. The spatial 
coordinate x is in units of \l h/muJT- 
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